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The Green's function formalism in Condensed Matter Physics is reviewed within the equation of 
motion approach. Composite operators and their Green's functions naturally appear as building 
blocks of generalized perturbative approaches and require fully self-consistent treatments in order to 
be properly handled. It is shown how to unambiguously set the representation of the Hilbert space 
by fixing both the unknown parameters, which appear in the linearized equations of motion and 
in the spectral weights of non-canonical operators, and the zero-frequency components of Green's 
functions in a way that algebra and symmetries are preserved. To illustrate this procedure some 
examples are given: the complete solution of the two-site Hubbard model, the evaluation of spin 
and charge correlators for a narrow-band Bloch system, the complete solution of the three-site 
Heisenberg model, and a study of the spin dynamics in the Double-Exchange model. 
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I. INTRODUCTION 

The physical system analyzed in this paper is an ag- 
gregate of interacting W^annier-electrons living on a lat- 
tice spanned by the vectors i. For the sake of simplic- 
ity, we restrict our study to single-band electron models; 
the generalization to multi-band models is straightfor- 
ward. The system is enclosed in a finite, but macroscop- 
ically large, volume V, containing M sites of the lattice, 
and is supposed to be in thermodynamic equilibrium at 
a temperature T. In a second-quantization scheme the 
dynamics of this system is ruled by a certain Hamilto- 
nian H = H [ip(i)] describing, in complete generality, the 
free propagation of the electrons and all the interactions 
among them and with external fields (e.g., electromag- 
netic fields, pressure and temperature gradients,...). cp(i) 
denotes an Heisenberg electronic field [i — (i, t)] in spino- 
rial notation satisfying canonical anticommutation rela- 
tions. Any physical property of this system can be con- 
nected to the expectation value of a specific operator. 

The expectation value (A) of any operator A = A [<p(i)] 
can be computed, for the grand-canonical ensemble, by 
taking the normalized trace of the operator weighted with 

the quantum-mechanical statistical factor . 
N = J2ia Vo-(*)Vo-(*) is the total number operator, (3 is 
the inverse temperature and fi is the chemical potential, 
which is fixed in order to get the desired average num- 
ber of particles N = (N). The chemical potential will 
be a function of N and T, as well as other parameters 
present in the Hamiltonian. Although the trace can be 
taken over any basis (i.e., over any complete set of states 
in the Hilbert space of the system) , the most convenient 
one, the eigenbasis, is constituted by the simultaneous 
eigenstates of H and N. If such a basis is known, then 
all the properties of the system can be exactly calculated: 
this procedure is known as exact diagonalization (ED). It 
is worth reminding that the Hilbert space of a fcrmionic 
system contains only those states compatible with the 
Pauli principle (i.e., states with occupation numbers per 
site and spini equal to either or 1). 

Generally, ED can be effectively applied only to sys- 



tems that are non-interacting or interacting, but very 
small. In particular, if the system is non-interacting the 
eigenbasis coincides with the canonical basis of the Fock 
space of the system (i.e., the set of states constructed 
by locating the electrons, one at a time, on the lattice 
sites in accordance with the Pauli principle). For small 
systems, it is always possible to exactly diagonalize the 
Hamiltonian according to the reasonable small number 
of available states, but when large interacting systems 
are considered the number of states can be enormous 
and ED is practically not applicable. This considera- 
tion gave birth to numerous numerical techniques: Lanc- 
zos, quantum Monte Carlo,..., which can be considered 
as attempts to construct an approximate version of ED 
that could be applied to very large systems. However, 
these numerical techniques have some very severe lim- 
itations coming from the unavoidable small number of 
sites they can treat (the computational time increases 
exponentially with the number of available states): they 
cannot give a reliable description of systems with long 
range interactions; phases presenting long range order of 
any kind are absolutely unaccessible; the very low res- 
olution in frequency and momentum prevents the ap- 
plications to systems with relevant low-energy features 
(e.g., systems that present Kondo-like effects) or with 
strong spatial-dependence or anisotropy in their physi- 
cal properties (e.g., systems that have a Fermi surface 
ill-defined, nodal or with high angular-momentum sym- 
metry). Moreover, the information we get by means of 
these techniques for a system of a certain size difficultly 
can be used for a system of bigger size and, even worse, 
does not give any clear idea of what can happen in the 
corresponding bulk system. 

According to this, we have to find an alternative exact 
analytical technique that can generate, for large interact- 
ing systems, approximate treatments not suffering from 
the very severe limitations noticed in the numerical meth- 
ods. In principle, this technique will obviously give the 
same exact results of ED. Coming back to our original 
problem, the evaluation of the expectation value (A), it 
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is possible to use the equation of motion 



(1.1) 



in order to derive one or more equations for this quantity 
or, better, for the corresponding Green's function (see 
next section). Actually, the equation of motion (|1.1|) 
naturally generates higher-order operators (i.e., op- 
erators constituted by more and more elementary fields, 
some of them centered on farther and farther sites from 
i). The process can be iterated by time-differentiating 
the newly generated operators and a chain of equations 
of motion can be constructed. The obtained system of 
equations of motion closes on a complete set of eigen- 
operators of the Hamiltonian 



V(M),tf 



e(i,j)V(j,t) (1-2) 



where ip(i) is a n-component spinorial field and e(i,j), 
usually called the energy matrix, is a square matrix of 
rank n. This approach is known as the equations of 
motion method (EM) and can be applied, obviously 
giving the same exact results, in all cases where we can 
also apply ED : if the system is non- interacting the origi- 
nal electronic operators ip(i) are the eigenoperators of the 
Hamiltonian; for small systems the number of equations 
of motion to be solved simultaneously, in order to find 
the complete set of eigenoperators, is reasonably small 
and makes the application feasible. When large interact- 
ing systems are considered the number of eigenoperators 
rapidly increases (diverges in the thermodynamic limit) 
and EM cannot be effectively applied just as ED could 
not be. However, the main difference between the two 
procedures is that EM can be still used in some approx- 
imation not subject to the severe limitations noticed in 
the numerical techniques derived from ED. 

Any approximation derived from EM is based on some 
of the peculiar properties of eigenoperators (some of 
them are reported below). These properties are obvi- 
ously not enjoined by eigenstates that have to be con- 
sidered always as a whole (symmetry considerations can 
only reduce a brute force diagonalization to a more re- 
fined block diagonalization which is in any way unfeasi- 
ble in the thermodynamic limit). The iterated process 
of time-differentiation generates more and more delocal- 
ized eigenoperators in direct space (i.e., eigenoperators 
containing original fields siting on more and more dis- 
tant sites), which are less and less relevant as they have 
eigenenergies rapidly decreasing with the spatial size 
of the eigenoperator (i.e., with the maximum distance 
among the sites where the constituting original fields are 
sited). Although the total number of eigenoperators is 
equivalent to the number of possible transitions among all 
the eigenstates and, therefore, goes as this latter number 
squared (i.e., if the number of eigenstates is n, then the 
number of eigenoperators is 2 )' ^° study a specific 
physical property we need only to analyze the dynam- 
ics of the few eigenoperators relevant to it. Furthermore, 



the eigenoperators can be easily generalized to any size of 
the system and the dynamics of all sites can be studied at 
once; this is impossible for the eigenstates. According to 
this, for very small clusters too, where the application of 
ED requires undoubtedly less effort than that of EM, EM 
solution is preferable as it has the fundamental property 
to be scalable (i.e., it gives a lot of information about 
both EM solution of bigger clusters and the approximate 
EM solution of the corresponding bulk system) . 

The line of thinking described so far follows the de- 
velopments of the condensed matter physics in the last 
decades. Both ED and EM try to diagonalize the Hamil- 
tonian under study, but in two different spaces. The for- 
mer searches for the eigenbasis within the Hilbert space 
of the system, the latter seeks an operatorial basis 
within the field space generated by the application of the 
Hamiltonian to the original field and to its bosonic ag- 
gregations (i.e., to fields constructed by an even number 
of original fields). While the states of a system drasti- 
cally change with its size (i.e., the corresponding Hilbert 
spaces do not overlap), the operators just increase in 
number and complexity (i.e., the new field space include 
the old one). Moreover, the relevance of an eigenopera- 
tor, which is measured by the magnitude of the scale of 
energy it describes, usually survives any change in size 
of the system. We can also define as minimal cluster 
the smallest one allowing all the terms of the Hamilto- 
nian to act properly. Only eigenoperators obtained for 
systems realized on clusters at least equal to the minimal 
one can be trusted and used to describe properties of the 
corresponding bulk systems. 

In order to construct any approximation scheme in the 
framework of EM, a convenient generalization of the con- 
cept of correlation function is provided by that of Green's 
function^ (GF). The latter has some advantages in the 
construction and solution of the equations that determine 
it. Moreover, the GF contain most of, practically all, 
the relevant information on the properties of the system: 
expectation values of observables, excitation spectrum, 
response to external perturbations, and so on. Differ- 
ent types of GF can be constructed; we will consider 
real-time thermodynamic GF where the thermal averag- 
ing process of the Heisenberg operators is performed over 
the grand-canonical ensemble. 

The traditional approximation schemes, often based 
on perturbative calculations, use as building blocks the 
non-interacting GF. The mean-field formulation, which 
corresponds to the linearization of the equation of mo- 
tion (ITU [i.e., ij^(i,t) = EjeCiJ^CM), where e(i,j) 
is now a scalar function], also belongs to this category. 
An intense study has been performed along this line and 
many techniques have been set up: perturbation expan- 
sions on the basis of Feynman diagrams, Dyson equation, 
Wick's theorem, and so on. It is worth noting that in or- 
der to describe phases with different symmetries, these 
schemes need to become self-consistent. 

All these techniques rely on the hypothesis that the 
interactions among the electrons are weak and can be 
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treated in the framework of some perturbation scheme. 
However, as many and many theoretical and experimen- 
tal studies have shown with more and more convincing 
evidence, all these methods are not adequate to treat 
strongly correlated electron systems (SCES) and differ- 
ent approaches must be considered. In these systems, the 
fundamental concept of the electron as a particle with 
some well-defined properties completely breaks down. 
The presence of the correlations modifies the properties 
of the electrons and, at the macroscopic level, new parti- 
cles are observed, with peculiar properties entirely deter- 
mined by the dynamics and the boundary conditions 
(i.e., all the elements characterizing the physical situa- 
tion we wish to study). These new objects appear as the 
final result of the modifications imposed on the electrons 
by the interactions and contain, by the very beginning, 
a relevant part of the effects of correlation. 

As simple, but significative, example, let us consider 
an atomic model with a local interaction U between the 
electrons (i.e., H — U(p^<pf(p^(pi). This model is exactly 
solvable in terms of the Hubbard operators 

t-[i-Vrf» 

v = |pVJ f 

Due to the presence of the local interaction U, the origi- 
nal electrons ip(i) are no more observables and new stable 
elementary excitations, described by the field operators 
£(«) and rj(i), appear. 

On the basis of this evidence, one can be induced to 
move the attention from the original fields to the new 
fields generated by the interactions. The operators de- 
scribing these excitations can be written in terms of the 
original ones and are known as composite operators. 
Several approaches have been formulated where compos- 
ite fields are used as operatorial basis for developing ap- 
proximation schemea 3 i 4 i 5 i 6 i 7 i 8 i 9 i 10 i 1:L i 12 i 13 i 14 . All these ap- 
proaches are very promising: some amount of the in- 
teraction is already present in the chosen basis and this 
permits to overcome the problem of finding an appro- 
priate expansion parameter. However, a price must be 
paid. In general, the composite fields are neither Fermi 
nor Bose operators, since they do not satisfy canonical 
(anti) commutation relations, and their properties, be- 
cause of the inherent definition, must be self-consistently 
determined. They can only be recognized as fermionic or 
bosonic according to the number of constituting original 
particles. 

New techniques of calculus have to be used in order to 
deal with composite fields. In developing approximation 
schemes where the building blocks are now the propaga- 
tors of composite operators, one cannot use the standard 
version of the consolidated schemes; diagrammatic ex- 
pansions, Wick's theorem and many other prescriptions 
are no more valid for composite operators. There have 
been attemptsi^*iS*ii to extend these schemes, but al- 
though very good results have been obtained for spin 
operators 1 ^, the complexity (and often the ambiguity) of 



the analytical calculations required by the Hubbard op- 
erators (the simpler among the fermionic composite oper- 
ators) does not allow, at least at the present, an effective 
application of such techniques to real problems. The for- 
mulation of the GF method must be revisited. As it will 
be shown below, three serious problems arise when we 
wish to study the propagators of composite fields: 

1. the appearance of some unknown parameters as 
correlation functions of field operators not belong- 
ing to the chosen operatorial basis; 

2. the appearance of some zero-frequency constants 
(ZFC) as a consequence of the existence of zero- 
frequency modes; 

3. the necessity of fixing the representation where the 
GF are formulated. 

In most of the approaches found in the literature the 
solution to the previous problems is the following. 

^ In order to determine the unknown parameters sev- 
eral methods (arbitrary ansatz, decoupling schemes, use 
of the equation of motion,...) have been considered in the 
context of different approaches (Hubbard I approxima- 
tion, Roth's method, projection method, spectral density 
approach,...). All these methods suffer from the severe 
limitation of not being fully self-consistent. On the other 
hand, any approach based on the correct use of compos- 
ite operators is, by construction, a fully self-consistent 
approach. As shown in Ref. fl3l in the context of the 
Hubbard model, all these procedures lead to a series of 
unpleasant results: several sum rules and the particle- 
hole symmetry are violated, there is no presence of a Mott 
transition, all local quantities strongly disagree with the 
results of the numerical simulation. 

13 Any symmetry enjoined by the Hamiltonian induces 
a degeneracy among the eigenstates of the system. The 
equivalence of two or more eigenenergies implies the pres- 
ence of zero-energy modes. In the case of bosonic Green's 
functions these modes give rise to some unknown quan- 
tities that we will call ZFC. The ZFC are really relevant 
quantities as they are connected to fundamental physi- 
cal properties such as the compressibility and the spe- 
cific heat: they can be considered as a measure of the 
fluctuations, quantum and/or thermal ones, present in 
the thermal averages of the generators of the symmetry 
group, which are usually bosonic. The ZFC are usually 
fixed by requiring the ergodicity of the dynamics of the 
relative operators with respect to the Hamiltonian under 
study. This is clearly a very strong assumption. As it 
will be shown in the third section of this paper, there are 
non-trivial examples of exactly-solvable systems where 
the ZFC do not assume their ergodic value: if we would 
have forced the ZFC to assume it, this would have im- 
plied a zero compressibility, specific heat,... Furthermore, 
although the response functions do not explicitly depend 
on them, there is an implicit dependence due to the in- 
herent self-consistency of the entire scheme. According 
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to this, in general, these quantities must be calculated 
case by case. 

121 The knowledge of the Hamiltonian and of the oper- 
atorial algebra is not sufficient to completely specify the 
GF. The GF refer to a specific representation (i.e., to 
a specific choice of the Hilbert space) and this informa- 
tion must be supplied to the equations of motion that 
alone are not sufficient to completely determine the GF. 
As well known, the same system can exist in different 
phases according to the external conditions; the existence 
of infinite incquivalent representations^ where the equa- 
tions of motions can be realized, allows us to pick up, 
among the many possible choices, the right Hilbert space 
appropriate to the physical situation under study. The 
construction of the Hilbert space where the GF are real- 
ized is not an easy task and is usually ignored. The use 
of composite operators leads to an enlargement of the 
Hilbert space by the inclusion of some unphysical states. 
As a consequence of this, it is difficult to satisfy a priori 
all the sum rules and, in general, the symmetry prop- 
erties enjoined by the system under study. In addition, 
since the representation where the operators are realized 
has to be dynamically determined, the method clearly 
requires a process of self-consistency. 

In the Composite Operator Methodic (COM), as il- 
lustrated in the next Section, the three problems are not 
considered separately but they are all connected in one 
self-consistent scheme. The main idea is that fixing the 
values of the unknown parameters and of the ZFC im- 
plies to put some constraints on the representation where 
the GF are realized. As the determination of this rep- 
resentation is not arbitrary, it is clear that there is no 
freedom in fixing these quantities. They must assume 
values compatible with the dynamics and with the right 
representation. Which is the right representation? This 
is a very hard question to answer. From the algebra it 
is possible to derive several relations among the opera- 
tors (e.g., iptr(i)(p<r(i) = 0): we will call them Algebra 
Constraint relations (AC). This set of relations, valid at 
microscopic level, must be satisfied also at macroscopic 
level (i.e., when the expectations values are considered; 
e.g., {(Pa(i)fa(i)) — 0). We also note that in general 
the Hamiltonian has some symmetry properties (e.g., ro- 
tational invariance in coordinate and spin space, phase 
invarian.ee, gauge invariance,...). These symmetries gen- 
erate a set of relations among the n-point Green's func- 
tions: the Ward-Takahashi relations^ (WT). It is worth 
noting that many approximations present in the liter- 
ature do not fulfill these consistency requirements and, 
consequently, obtain wrong results. Now, certainly the 
right representation must be the one where all relations 
among the operators satisfy the conservation laws present 
in the theory when expectation values are taken (i.e., 
where all the AC and WT are preserved). Then, we im- 
pose these conditions and obtain a set of self-consistent 
equations that will fix the unknown correlators, the ZFC 
and the right representation at the same time, avoid- 
ing the problem of uncontrolled and uncontrollable de- 



couplings, which affects many different approximation 
schemes and has been here definitely solved. This is the 
main ingredient of the COM, together with the recipes^ 
that we have developed in the last years in order to choose 
the appropriate operatorial basis of composite operators 
according to the specific system under analysis. As re- 
gards this last issue, we wish to drive the attention on 
the procedure we propose as it can be considered a sys- 
tematic attempt to seek and build up (exact as much as 
it is possible) operatorial basis for interacting systems. 
This is a new frontier in condensed matter theory (quan- 
tum Hall effect, heavy-fermion systems, quantum critical 
points, competing unconventional ordering phenomena, 
breakdown of Fermi liquid picture, connections among 
spin, charge, orbital and lattice degrees of freedom, ...) 
and our procedure should be regarded as an attempt to 
revisit the established picture for strongly correlated sys- 
tems. 

The second section is devoted to revisit the GF for- 
malism in presence of composite fields and to estab- 
lish the COM as a general procedure to compute GF 
of highly correlated systems. In the third section of 
the paper we will illustrate the formalism by consider- 
ing some specific examples: the two-site Hubbard model, 
the three-site Heisenberg model, a narrow-band Bloch 
system in presence of an external magnetic field and the 
double-exchange model. For the two-site Hubbard model 
we compute the fermionic GF independently from the 
bosonic one by means of the AC. The latter also allow 
us to fix the ZFC of the bosonic GF, which result in 
not being ergodic, and to get straightforwardly the right 
representation. The solution of the tree-site Heisenberg 
model shows the impossibility to get any spontaneously 
ordered state at finite temperature in a finite system as a 
consequence of internal consistency in the proposed for- 
mulation. Moreover, it is really relevant the existing re- 
lation between the number of ZFC appearing in the GF 
and the presence of the magnetic field. In the case of 
a narrow-band Bloch system in presence of an external 
magnetic field we will see that the ZFC relative to the to- 
tal number operator, which is an integral of motion, has 
a non-ergodic value, even if we have an ergodic charge 
dynamics. The double-exchange model finally gives us 
the possibility to show one way to apply the proposed 
formulation to large interacting systems. In this case, 
we also show how to recognize the manifestation of the 
Mermin- Wagner theorem^ within this formulation. For 
the exactly-solvable models used in the examples (i.e., the 
Hubbard and the Heisenberg models), we have checked, 
although it could not be otherwise, that the proposed for- 
mulation reproduces the exact results coming from ED. 
Actually, the COM has been developed just for systems 
large and interacting; the applications to systems that 
are small or non-interacting have only to be interpreted 
as mere demonstrations of all features of the method and 
of its power and correctness. Finally, in the Appendices, 
we give: a generalized perturbative approach for strongly 
interacting systems; part of the derivation of the general 
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formulation; the derivation of the zero-temperature for- 
mulation; the GF expressions, dispersion relations and 
sum rules where the presence of the ZFC is explicitly 
taken into account. 



II. GENERAL FORMALISM 

This section is devoted to revisit the GF formalism in 
presence of composite fields and to establish the COM 
as a general procedure to compute GF of complex in- 
teracting systems. Owing to the difficulties in dealing 
with composite operators, reported in detail in the pre- 
vious section, the study is performed completely within 
EM. We start by considering a set of composite fields, 
chosen according to a well-defined recipe 20 . The fields 
can be of fcrmionic or bosonic nature, according to the 
physical properties we wish to study 22 . In the case of 
fermionic operators it is intended that we use the spino- 
rial representation. The set ip(i) satisfies a linear system 
of equations of motion [see Eq. 11. 2|) ]. If the fields ip(i) 
are eigenoperators of the total Hamiltonian, the equa- 
tions of motion are exact. Several examples will be given 
in Sec. III. If the fields ip(i) are not eigenoperators of 
the Hamiltonian, the equations of motion are approxi- 
mated and all the formalism is developed with the aim 
of computing and using the propagators of these fields as 
a basis to set up a perturbative scheme of calculations. 
In Appendix A, we give a sketch of a generalized per- 
turbative approach based on a Dyson equation (Eq. IA5|) 
designed for formulations using composite fields. Then, 
the total weight of the self-energy corrections is bounded 
by the weight of the residual source operator 8J(i) [see 
Eq. (|A2|) ]. According to this, it can be made smaller and 
smaller by increasing the components of the basis 4>(i) 
(e.g., by including higher-order composite operators ap- 
pearing in 6J(i)). The result of such procedure will be the 
inclusion in the energy matrix of part of the self-energy 
as an expansion in terms of coupling constants multiplied 
by the weights of the newly included basis operators. In 
general, the enlargement of the basis leads to a new self- 
energy with a smaller total weight. However, it is neces- 
sary pointing out that this process can be quite cumber- 
some and the inclusion of fully momentum and frequency 
dependent self-energy corrections can be necessary to ef- 
fectively take into account low-energy and virtual pro- 
cesses. According to this, one can choose a reasonable 
number of components for the basic set and then use 
another approximation method to evaluate the residual 
dynamical corrections (e.g., specially adapted versions of 
the non-crossing approximation or of the FLEX). 

By considering two-time thermodynamic GF 2 - 23-24 , let 
us define the causal function 

-vefa-tdWwm) (2-i) 



the retarded and advanced functions 



G% A (i,j) = ±9[±(t i -t j )](fy(z),^(j)] ri 



and the correlation function 



(2.2) 



(2.3) 



Here rj = ±1; usually, it is convenient to take T) = 1 
(ri = —1) for a fermionic {bosonic) set tp(i) (i.e., for a 
composite field constituted of an odd (even) number of 
original fields <p(i)) in order to exploit the canonical anti- 
commutation relations satisfied by ip(i); but, in principle, 
both choices are possible. Accordingly, we define 

_ f{AB} = AB + BA for„=l 
l ' ir >~\[A, B]=AB-BA forry = -l [ ' 

< • • • > denotes the quantum-statistical average over the 
grand canonical ensemble. From the equation (|1.2|) for 
the set ip(i), the Fourier transforms of these functions sat- 
isfy the following equations (we consider a translational 
invariant system) 



[ W - £ (k)]4% >A (k,a;) = /W(k) 
[co-£(k)]C(k,uO = 

where 

/W(k)=^([^(i,t),v t a*)] 



(2.5a) 
(2.5b) 

(2.6) 



is known as the normalization matrix. T indicates the 
Fourier transform. The most general solution of equa- 
tions 1)2. 5fl is 



G 



(v) 

C.R.A 



1=1 



(k) 



uj - uji (k) _ 
i7ra[ W - Wl (k)] S g£ A (k)} (2.7a) 



C(k,w)=£*[w-w,(k)] c«(k) 



(2.7b) 



i=i 



Sc raO*-) an( i c^(k) are not fixed by the equations of 
motion and have to be determined by imposing the ap- 
propriate boundary conditions. w;(k) are the eigenvalues 
of the matrix e(k). (k) are the spectral density func- 
tions and can be expressed in terms of the matrices e(k) 
and (k) as 



<7^ ) (k) = ^(k)E^ 1 ( k )4?( k ) 



(2. 



where f2(k) is the n x n matrix whose columns are the 
eigenvectors of the matrix e(k). The summations run 
over the number of eigenvalues of e(k) and V represents 
the principal value. 
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By recalling the boundary conditions G^\t < 0) = C(k, w) = 2n S(ui)T(k) 



and (t > 0) = it is immediate to see that 
s f(k) = - s W'(k) = >')(k) 



(2.9) 



Then, the retarded and advanced GF are completely de- 
termined in terms of the matrices e(k) and /^(k). 

The determination of '^(k) and c^(k) require some 
more work. On the basis of the calculations reported 
in App. B, it is straightforward to obtain the following 
results 



£ r( 1 . , )(k) = (i + „)r(k) 

ze-4(k) 

2tt 



c »(k) 



1 liie-^*) 



(2.10a) 
(2.10b) 

<r^)(k) VZeB(k) (2.10c) 



^°(k) = ^z^- W3 (k) VZeS(k) (2.l0d) 

where «4(k) and B(k) are explicitly defined in Eq. 1B4|) 
and r(k), the zero- frequency function, is defined as 



r(k) = J- £ cWfk) 

We see that Eq. I|2.10b|> requires that 
£ ^ ( - 1 '°(k) = 



(2.11) 



(2.12) 



This condition comes from the requirement that the cor- 
relation function in direct space should not diverge: a so- 
lution with X!zg^4(k) o^ -1 '^ (k) 7^ implies a divergence of 
the Fourier coefficients (k) for any finite temperature. 
This is admissible only if the divergence is integrable and 
the corresponding direct space correlation function re- 
mains finite. A finite value of YlieACk) a (k) is gener- 
ally related to the presence of long-range order (i.e., sym- 
metry breaking) and the previous statement is nothing 
but the Mermin- Wagner theorem^!. A detailed analysis 
of this point will be illustrated in Section III by investi- 
gating the Heisenberg and Double Exchange Models. 

By putting Eqs. (|2~3j) and l|2~TDj l into Eqs. (|2~7|) we get 
the following general expressions for the GF 



(k) 



uj — u>i (k) ± i S 



(2.13a) 



G^(k >w ) = r(k) 



i 



v 



iee(k) 



r (v,l) 



(k) 



r\ e 



-0w 



uj + i 5 uj — i S 
1 



// e 



-/3 uj 



uj — uJi(k) + i5 



uj — uJi(k) — i8 

(2.13b) 



2ir £ *[w-w,(k)]- 
iee(k) 



cA^k) 



^ e -^o;((k) 



(2.13c) 



As shown in Appendix C, Eqs. (|2.13|l hold also in the 
limit of zero temperature (i.e., in the limit (3 — ► oo). From 
these expressions it is possible to get dispersion relations 
and sum rules that take explicitly into account the pres- 
ence of the zero-frequency function (see Appendix D). 

We see that the general structure of the GF is remark- 
ably different according to the statistics. For fermionic 
composite fields (i.e., when it is natural to choose r] = 1) 
the zero-frequency function T(k) contributes to the spec- 
tral function, it is directly related to the spectral density 
functions by means of equation l|2.10bjl and its calcula- 
tion does not require more information. For bosonic com- 
posite fields (i.e., when it is natural to choose r\ = — 1) the 
zero-frequency function does not contribute to the spec- 
tral function, but to the imaginary part of the causal GF. 
The causal and retarded (advanced) GF contain differ- 
ent information and the right procedure of calculation is 
controlled by the statistics. In particular, in the case of 
bosonic fields (i.e., for rj = — 1) one must start from the 
causal function and then use 



G<r2(k, w )' 



3? 



G^iKuj) 



± tanh 3 
2 



G^iKoj) 



(2.14) 



C(k,w) 



1 + tanh 



/3u 



G 



Kuj) 



On the contrary, for fermionic fields (i.e., for r\ = 1) the 
right procedure for computing the correlation function 
requires first the calculation of the retarded (advanced) 
function and then the use of relations identical to those 
of Eqs. 12.141 but with the subscript R, A and C inverted 
and the minus sign in the last equation changed to =F. 

Moreover, it is worth noting that T(k) is undetermined 
within the bosonic sector (i.e., for rj = —1) and should 
be computed in the fermionic sector (i.e., for 77 = 1) by 
means of equation (|2.10bl) or equivalcntly by means of 
the following relation 



r(k) 



1 



lim wG^M) 



2 w->o 



(2.15) 



However, the calculation of cr^'^k) requires the calcula- 
tion of /^(k) that, for bosonic fields, generates unknown 
momentum dependent correlation functions whose deter- 
mination can be very cumbersome as requires, at least in 
principle, the self-consistent solution of the integral equa- 
tions connecting them to the corresponding Green's func- 
tions. In practice, also for simple, but anyway composite, 
bosonic fields the T(k) remains undetermined and other 
methods rather than equation l|2.10b|l should be used. 
Similar methods, like the use of the relaxation function^, 
would lead to the same problem. 
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The zero-frequency function T(k) is known in the 
riteraturo 25 i 26 i 27 i 28 i 29 i 30 i 31 as an indicator of the ergodic 
nature of the dynamics of the operator with respect 
to the Hamiltonian H . We recall that a quantity A has 
an ergodic dynamics if and only if 



lim (A(t)A) = (AY 

t — »oo 



(2.16) 



that is, if and only if its auto-correlation attenuates in the 
time. We have not to forget that the condition H2.16|l is 
the same as the standard ergodic requirement (i.e., equiv- 
alence of averages taken in time and over the phase space) 
only for statistical averages computed in the microcanon- 
ical ensemble 25 ; in other ensembles it holds only in the 
thermodynamic limit. By recalling the general expres- 
sion l|2.13cf) for the correlation function, the condition of 
ergodic dynamics for tjj(i) is 



_L£ e ik(i-j )r(k)HVj(i)) ^t (j) ) 



(2.17) 



It is worth noting that T(k) generally does not assume its 
ergodic value [i.e., that required by Eq. H2.17(l ] and has 
to be computed case by case according to the dynamics 
and the boundary conditions. For instance, for any finite 
system the statistical ensembles are not equivalent and 
the criterion (|2.16|l holds only in the microcanonical one. 
Moreover, the condition (|2.16|l is not satisfied by any 
integral of motion or, more generally, by any operator 
that has a diagonal part with respect to the Hamiltonian 
under studyS 7 - (i.e., by any operator that has diagonal 
entries whenever written in the basis of the eigenstates 
of the Hamiltonian under study). This latter consid- 
eration clarifies why the ergodic nature of the dynam- 
ics of an operator mainly depends on the Hamiltonian 
it is subject to. It is really remarkable that the zero- 
frequency constants (ZFC), which are the values of 
the zero-frequency function T(k) over the momenta for 
which A(k) 0, are directly related to relevant measur- 
able quantities such as the compressibility, the specific 
heat, the magnetic susceptibility,... For instance, we re- 
call the formula that relates the compressibility to the 
total particle number fluctuations 



K -P N 2 



(N 2 ) - N 2 



(2.18) 



According to this, in the case of infinite systems too the 
correct determination of the ZFC cannot be considered 
as an irrelevant issue (e.g., Eq. (|2.18|l holds in the ther- 
modynamic limit too). In conclusion, Eq. (|2.17|l gener- 
ally cannot be used to compute the ZFC. In the next 
section, we provide some examples of violation of the 
condition (|2.16|) . It is necessary pointing out, in order to 
avoid any possible confusion to the reader, that we are 
using (here and in the examples presented in the next 
section) full operators and not fluctuation ones (i.e., we 
use operators not diminished of their average value, in 
contrast with what it is usually done for the bosonic ex- 
citations like spin, charge and pair). According to this, 



the ZFC can be different from zero (i.e., be equal to the 
squared average of the operator), and still indicate an 
ergodic dynamics for the operator. 

Summarizing, by means of EM and by using the 
boundary conditions relative to the original definitions 
of the various GF we have been able to derive explicit 
expressions for these latter [see Eqs. <|2.13ll ]. However, 
these expressions can only determine the functional de- 
pendence of the GF: their knowledge is not fully achieved 
yet. According to the (anti)commutation relations, the 
normalization matrix I^ 1 ' (k) usually contains some un- 
known functions that have to be self-consistently calcu- 
lated together with the ZFC (and the energy matrix e(k) 
if we use some approximation scheme). These functions 
are static correlation functions (correlators since now on) 
of operators not belonging to the chosen basis. In princi- 
ple, one could introduce a new set of composite fields and 
repeat all scheme of calculations in order to calculate the 
unknown correlators. However, the new set will possibly 
generate other unknown correlators and the entire pro- 
cess of self-consistency might become very cumbersome 
and, in most of the cases, not convergent. An alterna- 
tive scheme of calculation can be proposed. Fixing the 
values of the unknown parameters and of the ZFC im- 
plies to put some constraints on the representation where 
the GF are realized. As the determination of this rep- 
resentation is not arbitrary, it is clear that there is no 
freedom in fixing these quantities. They must assume 
values compatible with the dynamics and with the right 
representation. Now, certainly the right representation 
must be the one where all relations among the opera- 
tors are systematically conserved when the expectation 
values are taken (i.e., where all the AC and WT are sat- 
isfied). It is then clear that a shortcut in the procedure 
of self-consistency can be introduced. We can fix the 
representation by requiring that 



(2.19) 



where the l.h.s. is fixed by the AC, the WT and the 
boundary conditions compatible with the phase under 
investigation and in the r.h.s. the correlation function 
C(k, lo) is computed by means of Eq. I|2.13c|l . Equa- 
tions H2.19fl generate a set of self-consistent equations 
which determine the unknown parameters (i.e., ZFC 
and unknown correlators) and, consequently, the proper 
representationiSii 3 ^ 3 ^. It is worth noticing that the num- 
ber of constraints generated by Eqs. (|2.19|l can be differ- 
ent from the number of unknowns parameters. Generally, 
the coincidence of these two numbers signals that the cho- 
sen basic set gives a reasonable description of the dynam- 
ics contained in the truncated EM. Condition (|2.19|) can 
be considered as a generalization, to the case of composite 
fields, of the equation that, in the non-interacting case, 
fixes the way of counting the particles per site, according 
to the algebra, by determining the chemical potential. 
According to this, the unknown correlators, coming from 
the non-canonical (anti)commutation relations, have not 
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be seen like obstacles as many analytical techniques do, 
but like a possibility to fix the representation and sat- 
isfy all the symmetry relations. Any approximation not 
using them to do so will surely fail in reproducing the 
physics of the system under study. It is worth noting, 
and the examples of the next section will show how, that 
by means of Eqs. (|2.19[1 is often possible to close one 
sector (i.e., fermionic, spin, charge, pair, ...) at a time 
without resorting to the opening of all or many of them 
simultaneously. Obviously, this occurrence enormously 
facilitates the calculations. Finally, it is worth noting 
that the entire process of self-consistency [i.e., the use of 
Eqs. (|2.19l) ] will affect all the GF at the same time and, 
therefore, all the physical properties of the system. For 
instance, the linear response of the system to an exter- 
nal perturbation (susceptibility, conductivity, ...) is de- 
scribed by two-time retarded GF—. Although these type 
of GF do not explicitly depend on the ZFC, there is an 
implicit dependence through the internal self-consistent 
parameters, that is the unknown correlators. 

In this section, we have presented the general frame- 
work of the COM, which results to be a general method to 
deal with composite fields and, consequently, with com- 
plex correlated systems. In the next Section, we will illus- 
trate this calculation scheme by considering some specific 
examples. 



III. EXAMPLES 



where 








m- 


= [l-n(i)]c(i) 




(3.4a) 


T}(l) - 

i \ j 


= n(i) c(i) 




(3.4b) 








(3.4c) 


Vs{i) 




+ ^K tQ (*hM 


(3.4d) 



We define ip a (i) = . otij V'(j) and use the spinorial no- 
tation for the field operators. n M (i) = c^{i)a^c{i) is the 
charge (/j = 0) and spin (/i — 1,2,3) operator; greek 
(e.g., /i, v) and latin (e.g., a, 6, k) indices take integer 
values from to 3 and from 1 to 3, respectively; sum over 
repeated indices, if not explicitly otherwise stated, is un- 
derstood; ct m = (l,<x) and cr M = (— l,c?); a are the Pauli 
matrices. In momentum space the field tp(i) satisfies the 
equation of motion 

i J^-(M) = e(fcMM) (3.5) 

where the energy matrix e(k) has the expression 




-fi-2ta{k) -2ta{k) -2t -2t \ 

U-n 2t 2t j 

At -fi + 2ta(k) Ua(k) 

2t 2ta(k) U - n J 



(3.6) 



A. The two-site Hubbard model 

The two-site Hubbard model is described by the fol- 
lowing Hamiltonian 



where the summation range only over two sites at dis- 
tance a from each other and the rest of notation is 
standard—. The hopping matrix iy is defined by 



(3.2) 



where a(fc) = cos(fca) and k = 0, n/a. 

We now proceed to study the system by means of the 
equation of motion approach and the GF formalism 33 
described in Sec. [HI A complete set of fermionic eigen- 
operators of H is the following one 



Straightforward calculations, according to the scheme 
traced in Sec. [HI show that two correlators 

A = (C(i)e^)~(v a ^)v\i)) (3.7) 
P= J »/»(<)) - (ct(*)<u(<) [cjWcjw] ) (3.8) 

appear in the normalization matrix /(k) = 
.F({V>(M), ^(j,t)}). Then, the GF depend on 
three parameters: /i, A and p. The correlator A can be 
expressed in terms of the fermionic correlation function 
C(i,j) = (ip(i)ip '(J)); the chemical potential \i can be 
related to the particle density by means of the relation 
ri = 2 [1 — Cn(i, %) — C-2i(i, i)]. The parameter p cannot 
be calculated in the fermionic sector; it is expressed in 
terms of correlation functions of the bosonic fields n^ii) 
and C|(i)c|(i). According to this, the determination of 
the fermionic GF requires the parallel study of bosonic 
GF. 

After quite cumbersome calculations, it is possible to 
see24 that a complete set of bosonic eigenoperators of H 
in the spin-charge channel is given by 




( B<f\i) 



(3.3) flW(i) = 



(3.9) 
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where 










B^(i 


) = c f (i)o- M c(i) 






(3.10) 


B^(i 


)= C t(i)^«J a (i) 


- C t a (i) 


CT M C(0 


(3.11) 




) = d fi (i)-d«(i) 


+4(0 


-4 a (0 


(3.12) 


B*f\i 


) = d lt (i)-d°(t) 


-4w 


+4 a (0 


(3.13) 




) = UO)-fZ(i) 


-4(0 


+ 4 Q (0 


(3.14) 




) = /m(0-/£(0 


+ 4(0 


-4 a (0 


(3.15) 



with the definitions: 

d M (o = e f (o^»? a (o 

/ (0 = -^(0*7(0-^(0^(0 
+ ^(0*7(0^(0^(0 

/a(0 = £ f (0 m <(i) - \ie abc n b (i) n a c {i) 
The field B^(i) satisfies the equation of motion 

i^-B^>(k,t) = K(k)B^(k,t) 
at 

where the energy matrix n(k) has the expression 



(3.16) 

(3.17) 
(3.18) 

(3.19) 



/ 





-2f 











\ 


-4t [l-a(jfc)] 





u 






















u 


2f 













C7 








2t 










8/ 











V 











8t 





0/ 



(3.20) 



The energy spectra are given by 
wi(fc) = -24-^2 [1 - «(fc)] 



w 2 (7e) 
w 3 (fc) 
cji(k) 
u 5 (k) 
uj 6 (k) 

where 
Ju = 



2t^/2 [1 - a{k)\ 
-U - 4Ju 

U + AJu 



yV 2 + 64i 2 - U 



(3.21) 

(3.22) 
(3.23) 
(3.24) 
(3.25) 
(3.26) 



(3.27) 



Straightforward calculations according to the scheme 
given in Sec. |H] show that the correlation function has 
the expression 



1 6 

f3u n (k) 



k n—1 



1 + tanh : 



f M (k) (3.28) 



where 

(0) = /or n = 3,4,5,6 

(tt) = coth^^^) Vn 



(3.29a) 
(3.29b) 



Owing to the fact that zero-energy modes appear for n = 
1, 2 and k = [cfr. Eq. (|3~?T1) ]. ZFC appear in the 
correlation functions 



r (Al) (0) = i]T/ 



(o) 



(3.30) 



In principle r^(0) could be calculated by means of 
Eq. I)2.10b(l : however this would require the calculation 
of the anticommutators ({B^(i,t), B^ f (j,t)}) which 
generate correlation functions of higher order giving raise 
to a chain of GF whose closure is not evident. Similar 
methods, like the use of the relaxation function^, would 
lead to the same problem. One might think, as is often 
done in the literature, to fix this constant by its ergodic 
value. However, this is not correct as we are in a finite 
system in the grandcanonical ensemble and the ergodic- 
ity condition (|2.17|l does not hold. For the moment, we 
can state that this constant remains undetermined. 

The spectral density functions cr'™' M '(7c), calculated by 
means of Eq. (|2.8|) depends on a set of parameters which 
come from the calculation of the normalization matrix 
lM(k) = J 7 ([B<» (£,*), B^(j,t)]). In particular, for 
the (l,l)-component the following parameters appear: 



^ = (^(0^(0) 
c<* = (c a (i)J(i)) 

d = (cifyc^i) c\(i)c\(i) 
X a s = (n(i) • rT(z)) 
The parameters C a and Cf 2 



(3.31a) 
(3.31b) 

(3.31c) 
(3.31d) 

arc related to the fcrmionic 



correlation function C(i,j) = (ip(i) il>^(J)\ The parame- 
ter x" can be expressed in terms of the bosonic correla- 
tion function C^(i,j) = (B^(i) B^(j)). In order to 
use the standard procedure of self-consistency, we need 
to calculate the parameter d. For this purpose we should 
open both the pair channel and a double occupancy- 
charge channel (i.e., we will need the static correlation 
function (nt(i) n\ (i) n a (i))). The corresponding calcula- 
tions are reported in Ref. y3| where is shown that these 
two channels do not carry any new unknown ZFC. The 
self-consistence scheme closes; by considering the four 
channels (i.e., fermionic, spin-charge, pair and double 
occupancy-charge) we can set up a system of coupled self- 
consistent equations for all the parameters. However, the 
ZFC r^)(0) has not been determined yet: we have not 
definitely fixed the representation of the GF. 

In conclusion, the standard procedure of self- 
consistency is very involved and is not able to give a final 
answer because of the problem of fixing the ZFC. This 
problem is known in the literature as the zero-frequency 
ambiguity of the GF formalism 24 i 26 i 28 i 29 i 30 . 
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We will now approach the problem by taking a differ- 
ent point of view. The proper representation of the GF 
must satisfy the condition that all the microscopic laws, 
expressed as relations among operators must hold also 
at macroscopic level as relations among matrix elements. 
For instance, let us consider the fermionic channel. We 
have seen that there exists the parameter p, not explicitly 
related to the fermionic propagator, that can be deter- 
mined by opening other channels. However, we know 
that at the end of the calculations, if the representation 
is the right one, the parameter p must take a value such 
that the symmetries are conserved. By imposing the AC 
(|2.19l) and by recalling the expression for A we get three 
equations 



n = 2(1 - C u - C22) 
A = ~ C% 2 
C 12 = 



(3.32a) 
(3.32b) 
(3.32c) 



This set of coupled self-consistent equations will allow 
us to completely determine the fermionic GF. Calcula- 
tions show^ that this way of fixing the representation is 
the right one: all the symmetry relations are satisfied and 
all the results exactly agree with those obtained by means 
of ED. We do not have to open the bosonic channels; the 
fermionic one is self-contained. 

Next, let us consider the spin-charge GF. In the spin- 
charge sector we have the parameters C Q , Cf 2 , Xs > ^ an d 
the two ZFC 



^4E/n 0) (°) 

i—l 

&* = z£/i ( i fc) (°) fc = 1 > 2 > 3 



(3.33) 
(3.34) 




FIG. 1: bo and bk are plotted as functions of n for [7 = 4 and 
T = and 1. U and T are expressed in units of t. 
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Since we are in absence of an external applied magnetic 
field, bk takes the same values for any value of k. 

The parameter C a and Cf 2 are known, since the 
fermionic correlation functions have been computed. The 
parameters and d can be computed by means of the 
equations 



1 



»(*)) 



P 



X a s = (n(i)-n a (i)) 
The ZFC are fixed by the AC 

C[?{i,i) = { nti {i)nM 
By recalling (|3.28|) and (|3.29() we have 



1 6 



1 + coth 



with 



(?fy(i)7fy(i)) = 



n + 2D for = 
n-2D for fi = 1,2,3 



(3.35) 
(3.36) 

(3.37) 

on W 
(3.38) 

(3.39) 



FIG. 2: 60 and bk are plotted as functions of U for T = 0.01 
and n = 0.6, 0.8, and 0.9. U and T are expressed in units of 
t. 



D = (n^(i) ni(i)) is the double occupancy and can be cal- 
culated by means of the fermionic correlation functions 
D = n - 1 + C n . Eqs. lEOST and (pHgj) constitute a 
set of coupled self-consistent equations which will deter- 
mine completely the Green's function in the spin-charge 
channel. Calculations show that this way of fixing the 
representation is the right one: all the symmetry rela- 
tions are satisfied and all the results exactly agree with 
those obtained by means of ED. 

The ZFC bo and bk are plotted as functions of n and U 
in Figs, n an d 121 respectively, for various temperatures. 
It is worth noting that they assume their ergodic values 
(i.e. n 2 and 0, respectively) only in some regions of the 
parameter space: (at zero temperature) at n = 1 (both 
60 and bk) and at n = 0.5 (60 on ly)- I n these regions, 
the grand-canonical ensemble is equivalent to the micro- 
canonical one and the underlying ergodicity of the charge 
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and spin dynamics emerges. 

It is worth noting that the ZFC bo is directly related to 
the compressibility by means of the following relation 33 



1 



k^T n 2 



\b - n 2 



(3.40) 



According to this, if we erroneously set the value of bo 
to the ergodic one (i.e., n 2 ) we would get a constant zero 
compressibility. 



B. The three-site Heisenberg model 

The three-site Heisenberg model, in presence of an 
external magnetic field h, is described by the following 
Hamiltonian 



H = J £(•) • S a (*) - h S, (t 



(3.41) 



i=l 



i=l 



where S(i) is the local spin at the site i, with quantum 
number S = | • The relative positions of the three sites 
are those of the end-points and the middle-point of a seg- 
ment, in this case we are using periodic boundary condi- 
tions, or those of the vertices of an equilateral triangle; 
in both cases, the distances of two of them in the given 
order (i.e., 1 — > 2, 2 — > 3, 3 — > 1) is taken to be unitary. 
The notation S a (i) indicates 



S ?tt (i,<) = ^o ij 5(j,t) 
j 

The projection operator ay is defined by 

a ^E cik(i ~ j)a ( k ) 



a(k) = cos(k) 
where k = - 2 f, 0, \f 

' 'da 

A complete set of eigenoperators of H is 

V; (m) (z) 



(3.42) 

(3.43) 
(3.44) 



(3.45) 



where 



4 m) (*) = - 


r s+(i) = s x - 


f iS y for to 
for to 


= 1 
= 2 


(3.46) 


v4 m) (*) = - 


r = i x + 


ily for m = 
for m = 


1 

2 


(3.47) 


</4 m) (*) = - 


u + (i) — u x - 

L "*(*) 


- iu y for rn 
for to 


= 1 

= 2 


(3.48) 



Hereafter, we will use the two sets of indices {x, y, z} and 
{1,2,3} interchangeably. 



The composite fields lk(i) and Uk(i) are defined as 

lk(i) = ie kpq S2(i)S q (i) (3.49) 
u k {i) = ie kpq [l%(i)S q (i) + S$(i)l q (i)] (3.50) 

The field ^^{i) satisfies the equation of motion 

i4v> (m) (0 =£ (m V (m) W (3-51) 

where the energy matrix e^™) has the expression 

a TO /i 2J 
e ( m ) =| a m h 2J | (3.52) 
|J a m /i 



with a m = 1 — (52 m . The energy spectra are given by 

— a m h (3.53) 
uj^ m) = l(2a m h-3J) (3.54) 



aj^" 1 ^ = —(2a m h + 3 J) 



(3.55) 



By means of the equation of motion (|3.51|) , the corre- 
lation function 



CW(i.i) = (^ m) (* ,(m) tl'); 

= lY^- /^e ik(i - j) - iw( * l - tj) C (m) (k,w) (3.56) 
3 k 2lT J 

has the expression 

3 

C (m) (k, u) = ^2s[u- uj { n m) (k)]c ( "' m) (k) (3.57) 

71=1 

where the matrices c^ n,m ^(k) have to be calculated. 

Straightforward calculations according to the scheme 
given in Section II show that the correlation function is 
given by 



c (1) (u)4EE eik(i 



ik(i-j)-i W W(fe)(t 4 -t 3 ) 



n=l k 



x [l + coth( ]]^ (M) (k) (3.58) 



c«(i,j) = 6 w (ij) + jEE e ' k( *" J 



-j)-i"i 2, W(t-t J ) 



n=2 k 



x [l+cothf ^)]a("- 2 )(k) (3.59) 



where the zero-frequency function 

& (2) (u) = ^5E eik(i_j)c(1 ' 2) (k) 



(3.60) 
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appears owing to the presence of the zero-energy mode 
= 0. 

The spectral density functions cr^"' m ^(k), calculated by 
means of Eq. (|2.8|1 , have the following expressions 





(k) 


= A (1 ' m) (k)A (1) 




(3.61) 


a (2, m 


(k) 


= \^ m \k)A^ 




(3.62) 


a (3,™ 


(k) 






(3.63) 


^(l,m 


(k) 


= /(r ) (k)-y4" 


>(k) 


(3.64) 


^(2 : m 


(k) 


= 3/ 1 ( ' 2 " ) (k)-44™ 


(k) 


(3.65) 


\(3,m 


(k) 


= 3/(™ ) (k)+4/(™ 


(k) 


(3.66) 



These parameters are expressed in terms of the cor- 
relation function C^ m ^(i,j) and self-consistent equations 
are easily written by means of Eqs. H3.58J1 and Ij3.59|l . 
However, in order to close the set of equations we need 
to know the zero-frequency constant 



6 (2)a 
u ll 



i^ a (k)6^(k) = ^i^ a (k) c (V 2) (k) 



(3.80) 



This quantity undetermined within the bosonic sector, 
can be obtained, as proposed in Section II, by fixing the 
representation of the GF by means of Eq. I|2.19l) . In 
particular, the AC requires that 



with 
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The normalization matrix /( m )(k) 

J 7 ([V' (m) (i^),V' t(m) (j^)]): that haS tllC form 



/W(k) 



J}i'(k) 7^(k) J^(k) 

J#(k) J$(k) &!&>(k) 

^22 (k) (k) (k) 

r( 2 )/ 



lif(k) 







^if(k) 
I#(k) 



Ir( 2 ) 
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with 



(3.67) 
(3.68) 
(3.69) 



(3.70) 



(3.71) 



#00 = 


2M 




(3.72) 


#(k) = 


-[1 


-a(k)] (c[\ )a + 2C{? a ) 


(3.73) 


4^(k) = 


-[1 




(3.74) 


4 2 '(k) = 


-[1 


- a(k)] Q - |c&>» + id?- 





(3.75) 



(S+(i)l-(i)) = -M + (S°(i)S z (i)) + - (S +a (i)S-(i)) 

(3.81) 

This equation, together with the others coming from 
the definitions l|3.58|l and (|3.59() . gives a set of five cou- 
pled self-consistent equation for the five parameters M, 
C[^ a , C[^ a , Cj^, &n a ■ The system can be analytically 
solved. In particular, the magnetization per site M and 

(2) 

the zero- frequency constant b\i are given by 

+ 3 cosh(/?/i) 



1 , / P h 
M = - tanh — 

6 V 2 
( 2) Q _ 9 cosh(/3fe) - 



cosh(/%) 



36 e^' 



cosh(/3/i) 



(3.82) 
(3.83) 



It should be noted that other ZFC appear into the 
model. Again, they can be fixed by the Algebra Con- 
straint. For example, the ZFC 



»i 2 HE^( k )= 2^£^ 2) (k) 

k k 

is determined by means of the equation 



and takes the value 

(2) _ 9cosh(/3/i) -4 + 5e' 3 i 
"li — 



(3.84) 



(3.85) 



36 e^ J +cosh(/3/i) 



(3. 



depends on the set of parameters 



1 



M = (S z (i)) = -- = (S+(i)S-(i)) - 

c$ a = (s + «(i)s-(i)) 

C[f a = (S»(i)S z (i)) 



a 



(i) 

22 



(3.76) 

(3.77) 
(3.78) 
(3.79) 



We note that the two ZFC &xi and assume the 
ergodic value M 2 only in the limit of very large external 
magnetic field h. In Figs. 01 and 01 the zero- frequency 
constants o^ a and are plotted as a function of the 
magnetic field h/J and temperature T/J, respectively. 
For comparison the ergodic value M 2 is also reported. 

Let us now consider the S -channel (i.e., m = 1) and 
the relative ZFC. For h ^ \J there are no zero-energy 
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FIG. 3: The zero-frequency constants b^ a and are plot- 
ted as a function of the magnetic field for T/J — 0.1. For 
comparison, the ergodic value M 2 is also given. 
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FIG. 4: The zero-frequency constant fc^ and b 2 ^ are plotted 
as a function of the temperature for h/J = 1. For comparison, 
the ergodic value M 2 is also given. 



modes and the ZFC relative to the S ,± -operators assume 
the ergodic value, i.e., zero. For the special case h = | J 

a zero-energy mode uj^ — appears and the ZFC be- 
come nonergodic. Again, the Algebra Constraint can be 
used to fix these quantities. Straightforward calculations, 
according to the proposed scheme, give 



(k) 

3e^ J + 4 



l 8e -3/3J (i + ffiVy 



3e^l J + 1 



9e -3/3J (l + e /3|J^ 

Let us now consider the limit of zero temperature. 
From the previous expressions it is easy to derive the 
following. 



• Case J > (Antiferromagnetic exchange) 
M= {\ for^>|J 

[ \ for h<y 



The two ZFC b^ a and bf} take the values 



(3.87) 



6 (2) 

6 (2) 
°11 



t(2)a 



1 

4 



_ 5 , (2)a _ _5_ 

5V "36 



for h> - J 
2 

for h <-J 
2 



(3.88) 
(3.89) 



They are ergodic for h > | J and nonergodic for h < 

§ J . The other two ZFC b^ a and b$ take the ergodic 
value (0) for any h > 0. In the limit of h — > we obtain 
the ferromagnetic solution. 

• Case J < (Ferromagnetic exchange) 



(3.90) 



All the zero-frequency constants take the ergodic value. 
In the limit of h — > we obtain again the ferromagnetic 
solution. 

Let us now consider the case of absence of external 
magnetic field (h = 0). For this situation there are two 
energy modes, one in the 5* z -channel (i.e., m = 2) = 

and one in the S' ± -channel (i.e., m — 1) — 0. In 
order to avoid divergencies in the correlation functions, 
it must be cr^^^k) = for m = 1, 2 and for all values 
of k. It must be 



= 



T (m) 



(k) 







M = 
1 



(3.91) 



C$ a + 



2C[f c 



8C£> = 0(3.92) 



By solving the self-consistent equations and by means 
of the Algebra Constraint, one finds the following noner- 
godic values for the ZFC : 



A2) 
U ll 

fc (2)a 
°11 



-b W - - 



5_ 
3~6 



I 6 (l)a _ 



e^2 



36 (e' 3 !' 7 + l) 



(3.93) 
(3.94) 
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b {2)a 



Summarizing, we have 



for J > 
for J < 



(3.95) 
(3.96) 
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71/ 




S ± 
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E 
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1 
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N 
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1 
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E 





< l J 
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1 
fi 
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E 





>fj 
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1 
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E 












N 


N 


7^0 


—> oo 


7^0 


1 
2 


E 


E 


7^0 


^0,| J 


7^0 


7^0 


N 


E 


7^0 


^ 7 


7^0 


7^0 


N 


N 



where E and N stays for an ergodic and nonergodic be- 
havior of the corresponding ZFC, respectively. The first 
two lines of the table consider the cases in which we have 
a ferromagnetic solution. 

The operator S z — Xa^zW: as an y constant of mo- 
tion, has essentially a non ergodic dynamics. Actually, 
for high values of hj J, S z is forced to assume the higher 
possible value | and no fluctuations are allowed (i.e., the 
susceptibility vanishes): the dynamics returns to be er- 
godic. At zero temperature and for the ferromagnetic 
case (i.e., J < 0) the system is polarized and S z is er- 
godic for all finite value of the magnetic field and in the 
ferromagnetic phase. The operator S + = ^2iS + (i) nas 
an ergodic dynamics only in presence of the magnetic 
field h or in the ferromagnetic phase as it is no longer 
an integral of motion in these cases. Also, for the special 
case h = | J the operator S + becomes nonergodic. These 
results show how the ergodicity of the dynamics of an op- 
erator can strongly depend on the boundary conditions. 

It is worth noticing that the ferromagnetic phase has 
been obtained only at exactly zero temperature (i.e., 
when the applied magnetic field has been sent to zero 
after setting the temperature to zero). This is due to 
the size of the system; finite systems can sustain ordered 
phases only at exactly zero temperature. The correla- 
tion functions in direct space should be computed by fi- 
nite sums over momenta (see Eq. I3.56|) and for vanishing 
spectra (e.g., for vanishing applied magnetic field; see 
Eq. 13.53(1 the Bose factor (see the coth in Eq. 13.58(1 di- 
verges except at exactly zero temperature. Only in this 
latter case (i.e., T = 0) the corresponding spectral den- 
sity function can retain a finite value (and consequently 
the magnetization too; see Eas. 13.611 13.641 and 13.72(1 in- 
stead of being forced to vanish in order to avoid diver- 
gences in the direct space correlation functions. In prac- 
tice, we allow the magnetization to be finite and search 



for a fully self-consistent solution. The system will self- 
adjust by selecting only those states with a finite mag- 
netization of the same sign of that assigned as initial 
condition according to the ergodicity breaking inherent 
to any symmetry breaking. 

We wish to remark that all the results obtained in this 
section exactly agree with those obtained by means of 
ED. 



C. A narrow-band Bloch system in presence of an 
external magnetic field 

A narrow-band Bloch system in presence of an external 
magnetic field is described by the following Hamiltonian 



(3.97) 



where 713(1) is the third component of the spin density 
operator and h is the intensity of the external magnetic 
field. The indices i and j run on an infinite d-dimcnsional 
lattice. Straightforward calculations show that the causal 
Green's function G£'(i,j) — (T [n ^(i) n ^(j)]) and the 
correlation function C^(i,j) — (n^(i) n^(j)) of the 
charge-spin operator n^(i) = c* (i) c(i) have the fol- 
lowing expressions 

(k, u) = -i (2i:) d+1 a- d S (d) (k) S(u) r&> - Q (p) (k, u) 



(k, u) = (2ir) d+1 a- d S (d) (k) 5{w) 

Q^(k,«) 



1 + tanh — 



(3.98) 



(3.99) 



where 8^ d '(k) is the ci-dimensional Dirac delta function. 

comes from the proper fermionic loop and is 
the Fourier transform of 



(hj) = Tr [dp G c (i,j) o> G c (j, i)] 



(3.100) 



Here Gc(i,j) = (T [c(i) c^(j)] ) is the causal fermionic 
function and has the expression 



G, 



0E n (k.) 



LU-E n (k)+i5 uj-E n (k)-iS 



with 



£i(k) 
E 2 (k) 



-fi-2dta(k) - h 
-fi - 2dta(k) + h 





I 


: 




: 

















(3.101) 

(3.102) 
(3.103) 

(3.104) 
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where 



1 d 

a(k) = — cos(fci a) 



(3.105) 



The ZFC is fixed by the AC fFTfy which requires 



(27T) 



d+l 



cfjfcdw 



1 + tanh 



/3w 



^5 



(3.106) 



The loop QO)(k,w) can be calculated by means of 
(|3.101|l . Calculations show 



1 + tanh 



/3w 



{2TT) d + 1 

= ( n ) - ( n l) 2 - ( n l) 2 f° r M = 0, 3 
= <n)-2<n T (*))(ni(i)) for » = 1,2 



9[Q M (k,w)] 



(3.107) 
(3.108) 



By recalling the AC (JSSJl, Eq. (|3~TU^ gives for the 



rw - <n} 2 

r (l,2) =Q 

r< 3 > = (n 3 f 



(3.109) 
(3.110) 
(3.111) 



in accordance with the ergodic nature of the spin and 
charge dynamics in this system. 

It is worth noting that the compressibility of this sys- 
tem can be computed by means of the general formula 
(|2.18|) that holds in the thermodynamic limit too and 
gives 



1 a d 



<»> a 2 2 ^ d ti 



(k) 



(3.112) 



where C n (k) = cosh 2 ^ E ^^ ^j. We can see that an 

ergodic charge dynamics can lead to a non-ergodic value 
of the ZFC relative to the total number operator, which 
is an integral of motion. Also in the infinite systems 
the decoupling inspired by the requirement of ergodicity 
cannot always be applied. 

D. The Double Exchange Model 

The Double Exchange Model is defined by the follow- 
ing HamiltonianM 

H = E (*« - ct w Jh E ( 3 - 113 ) 

ij i 

s(i) is the spin density operator of the electron and is 
given by s(i) = ^c 1 (i) a c(i); S(i) is a localized spin; Jh 
is the ferromagnetic Hund coupling (Jh > 0). In the 



nearest-neighbor approximation for a d-dimensional cu- 
bic lattice with lattice constant a, iy takes the form 

« = - 2rf 4E eik ' (H>a(k) (3 - 114) 



tjj = —2dta\ 



where a(k) has been defined in the previous section and 
d is the dimensionality of the system. Let us introduce 
the Heisenberg field 



B(i) = ( S+ W 



(3.115) 



where s (i) and S (i) are the standard rising and low- 
ering spin operators. 

This field satisfy the equation of motion 



yJ dt y! l-J H A(i) 



where 

p(i) = c\ a (i)c l (i)-c\(i)c'{(i) 
\{t) = s+{t)S z {t)-s z {t)S+{t) 



(3.116) 



(3.117) 
(3.118) 



We linearize the equation of motion 13.1161) by projecting 
the source J(i) on the basis (|3. 115(1 



j(z)«5> B (ij)i?a<) 



(3.119) 



where the coefficients are determined by the following 
equation 

< [ j(i, t), St a «)] > = x; !) ( [^a. *)> Bt a. *)] ) 
i 

(3.120) 

Let us compute, within the framework described in 
Sec. |nj the causal Green's function 



G(k,u)=F(T [B(i)Bt{j)]) 
cr(»)(k) 



E 



p 



LO - w n (k) 



-i7r*[w-w„(k)]5W(k) 



(3.121) 



and the correlation function 



C(k, w) = ^ <B(») 5+ (i)> = X * [w - w„ (fc)] c<"> (k) 

(3.122) 



n=l 



5 (n) (k) and c^( k ) 

are still unknown functions; w n (k) 
are the eigenvalues of the matrix e B (k); crW(k) are the 
density spectral functions, completely determined by the 
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matrices e B (k) and I B {k) = T {[B(i,t), (j,t)]) by 
means of relation l|2.8|) . We have 



I B (k) 



2 («,(»)> 
2 <£,(<)) 



(3.123) 



For the sake of brevity, the explicit expressions for the 
energy matrix £ S (k), the energy spectra w„(k) and the 
spectral density functions c^(k) are not reported here 
and can be found in Ref. The calculations show that 



lim wo(k) = 



lim cr (2) (k) 



-'ll i 22 



k o /, ;> ; 



J 22 



1 -1 
-1 1 



(3.124) 
(3.125) 



According to the scheme of calculation given in Sec.lTTl 
we generally have 



cW(k) = 2 7 r<5 n2 <5(k)r i 



1 + coth 



/M k ) 



cj(")(k) 
(3.126) 

S< n >(k) = 25 n2 S(k)T B +coth^^a(")(k) (3.127) 
The Green's functions have the following expressions 
G(k,w) = -2i<J(o;)5(k) T B 



o-w(k) 



1 



w — Lu n (k) + iS lu — w n (k) — i <5 



The spectral density function cr( 2 )(k) vanishes, in 
agreement with the general relation 12.12fl . 



T = In this case Eq. (|3.128|l becomes 

(4 rB =< B(i)B,<i) > 

-jit frt«h(k)k l "(k) (3.131) 

where • ■] is the ordinary step function. We see that 
we can have ferromagnetic order for any dimension. 

The results of the calculations for this model illustrate 
what stated in Section II regarding broken symmetry 
states in bulk systems. When a zero-energy mode ap- 
pears in the ladder operator sector (e.g., S + ) the cor- 
responding spectral density function, which is related to 
the order parameter (e.g., a cx / cx (S z )), must vanish for 
all finite systems at finite temperatures and for all infinite 
systems when the divergence of the Fourier coefficients of 
the correlation function is not integrable in order to avoid 
divergences in the direct space correlation functions. At 
zero temperature it is always possible to have a finite 
value for the spectral density function and consequently 
for the order parameter. This result is a manifestation of 
the Mermin- Wagner theorem™ which prevents the sys- 
tem from approaching, in some dimension, a particular 
ordered phase except at zero temperature. 



C(k,w) = 2nS{uj)S(k)T I 

2 



71=1 



1 + coth 



/9w n (k) 



a (n) (k) 



The ZFC T B is determined by means of the A C which 
requires 



(2n) 



tT 1 



(B(i)Bt(i)) 



2 
\ ^ 



2(2*)" ^ 



d d k 



1 + coth 



0w«(k) 



An) 



(k) 



(3.128) 



In the case of a three-dimensional system the integral 
in Eq. (|3.128J) is finite. The Green's functions are fully 
determined and a ferromagnetic order does exist. In the 
case d < 3 we must distinguish two cases. 

T > In this case the divergence of the integrand in 
Eq. (|3.128l) is not integrable and the integral is di- 
vergent. The only physical solution is absence of 
ferromagnetic order. The magnetization must van- 
ish 



Zf 1 = 2( Sz (z)) = 



i£ = 2 (&(»)) 







(3.129) 
(3.130) 



IV. CONCLUSIONS 

In conclusion, the GF formalism for composite oper- 
ators has been revised by making use of the equations 
of motion method. It has been shown that all the gen- 
eral relations (spectral representation, sum rules, etc.) 
can be derived without resorting to the knowledge of the 
complete set of eigenstates of the Hamiltonian. The ad- 
vantage of using the equations of motion formalism is 
that it can be applied to any operatorial basis both exact 
and approximate. Special attention has been paid to the 
presence of the ZFC and to the problem of determining 
unknown parameters related to higher order correlators. 
The ZFC issue is quite relevant because such quantities 
are directly related to many response functions. We have 
shown that an effective and proper way to fix the repre- 
sentation is to impose the constraints coming from the 
AC and the WT. When these conditions are required, a 
set of self-consistent equations is obtained that permits to 
compute both the parameters appearing in the spectral 
functions and the zero-frequency component of the GF, 
avoiding the problem of uncontrolled and uncontrollable 
decouplings, which affects many different approximation 
schemes and has been here definitely solved. 

Moreover, it is worth reminding the following issues, 
which have been discussed in detail all over the text: 
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• The two-time retarded (advanced) and causal 
bosonic GF carry substantially different informa- 
tion. 

• The ergodicity condition cannot be used a priori to 
compute the ZFC. 

• The Mermin- Wagner theorem^ naturally appears 
as a requirement to avoid divergences in the direct 
space correlation functions. 

It is also necessary pointing out - we already did it in 
Section II - that, although a careful choice of the com- 
ponents for the basic set makes possible the description 
of the main scales of energy present in the system under 
analysis, the inclusion of fully momentum and frequency 
dependent self-energy corrections can be sometime nec- 
essary to take into account low-energy and virtual pro- 
cesses. 

The calculation scheme has been illustrated by con- 
sidering four systems: the two-site Hubbard model, the 
three site Heisenberg system, the narrow-band Bloch sys- 
tem and the Double-Exchange model. These examples 
clearly show the relevance and complexness of the above 
issues and illustrate in detail the application of the pro- 
posed procedure. It has been checked that the proposed 
scheme gives the exact result when solvable systems are 
considered. 
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APPENDIX A: GENERALIZED PERTURB ATIVE 
APPROACH FOR SCES 

Given a certain Hamiltonian H = H [(p(i)\, where ip(i) 
denotes an Heisenberg electronic field [i — (i, t)] in spino- 
rial notation satisfying canonical anticommutation rela- 
tions, and a set of composite operators chosen in the 
spirit of the discussion given in Section I, the equations 
of motion for the propagator of the field ip(i) can be ob- 
tained by the dynamics obeyed by this latter which reads 
as 



(Al) 



In complete generality, this equation can be rewritten 



as 



^(i) = EW) + ^) 



(A2) 



where the linear term e ip represents the projection of the 
source J(i) on the basis tp(i). The energy matrix e(i, j) 
can be computed by means of the equation 



i [47(i,t),V t a*)],)=0 (A3) 
which defines the residual source SJ (i) and gives 



j) = ]T ( [J(i, f),^(l, t)] n )( t), ^(j, t)] 



(A4) 

Obviously, also less systematic projections of the 
source could be attempted and will result in different 
determinations of e(i, j) and 6J(i). 

After Eq. (|X2l) , the Fourier transform Gq ) (k, w) of the 

GF Gg\i,j), where Q = R (retarded), A (advanced), 
C (causal) (see definitions in Section II), satisfies the fol- 
lowing equation 



G$> (k >W ) 



7 (?)) (k) 



E^k^G™^) (A5) 



where the propagator Gg (k, u>) is defined by the equa- 



tion 



[ W -e(k)] Gg> (k, W )=jM(k) 



(A6) 



The matrix iw(k), known as the normalization ma- 
trix, is defined as 



/(")(k)=^([^(i,t),V t (j,t)] 1) ) (A7) 
Eg ' (k, uj) is the proper self-energy and has the expression 



E^(k,c)=<i r (k, W )j(")(k)- 1 



(A8) 



where Bq \ rr (k, u) is the irreducible part of the propaga- 
tor B ( Q>\k,uj) =F(Q [5J(i)&P(j)]). Equation ESJ can 
be formally solved to give 



Gg>(k,a,) 



w-e(k) -E^k.w 



(A9) 



Equations (|A5|1 and (|A9(I are nothing else than the Dyson 
equation for a formulation based on composite fields and 
represents the starting point for a perturbative calcula- 
tion in terms of the propagator Gq^ (\s.,u>). The proper- 
ties and the determination of this latter are derived and 
discussed in Section II. 
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APPENDIX B: KMS RELATION AND THE 
GENERAL FORMULATION 

From the definitions H2.1(l - (|2.3|) we can derive the fol- 
lowing exact relations 



[m,^U)].r,) (Bla) 



%Hi,j)-G%\i,j) = ^(J)] v ) (Bib) 



G 



By making use of the Kubo-Martin-Schwinger (KMS) 
relation (A(t) B(t')) = (B(t') A(t + where A{t) 

and B(t) are Heisenberg operators at time t, the rj- 
commutator can be expressed in terms of the correlation 
function as 



[l + rje-^] C(k,w) (B2) 



where M is the number of sites and k runs over the first 
Brillouin zone. Then, the equations (|B1|) in momentum 
space become 



n 

X;*[ W - WI (k)] {<$ ,0 (k) 



~ [l -rye"^] c (i) (k)j = (B3a) 



f^^-^(k)]{^(k) 

1=1 

~ [l + Ve'^] cW(k)| =0 (B3b) 

In order to solve these equations, we have to take into 
account that for any given momentum k we can always 
write 

|=0 fi»iE^k)C« = {l n} 

w for / G B(k) = N - A(k) 

Obviously, .A(k) can also be the empty set (i.e., A(k) = 
and B(k) = N). Combined use of Eqs. lfB3|) and Ip4| 
gives the results reported in Sec. II [Eq. I|2.10[l ]. 



APPENDIX C: FORMULATION IN THE LIMIT 
OF ZERO TEMPERATURE 



At zero temperature Eq. (|B3(1 is not applicable and we 
should proceed in the following way (the usual derivation 



in terms of a complete set of eigenstates of the Hamilto- 
nian can be found in Ref. |3r|). 

Let us consider the correlation functions 



1>-^1 



N ^ 2tt 

A; 



d^jMi-i)-" (u-W] c t (k )W ) 



= ^Ei/^ ei[k - ( " j) ^^ ,] ^(k, W ) 

k 

By taking the limit T — ► of the KMS relation 

C^(k,Gj)=e-P"C^(k,Gj) (CI) 

it is immediate to see that for any finite value of the 
Fourier coefficients, it must be 



CWt (k, uj) 
C^^(k,oj) 



^ for uj > 
= for uj < 



1 = for uj > 
\ 7^0 for uj < 

Furthermore 

C^(k,0) = C^ t (k,0) 



(C2) 
(C3) 

(C4) 



Let us consider the energy spectra wj(k) and let us 
write in complete generality, for any given momentum k 



= for I G A(k) C H 
uj t (k) = <( > for I G C(k) C H 
< for I G V(k) C H 



(C5) 



Then, Eqs. (|B1|I in momentum space are written as 



«5M 53 U^w-^a-^cS^k)} 

ie.A(k) 1 '' J 

+ £ ^^-^(k)]{^ i3 (k)-^-c« t (k)} 
zeC(k) L ^ J 

+ J2 S l" ~ {-9c '°(k) + |;<$^k)} = 



iex>(k) 



(C6a) 



53 5[u; - WI (k)] ^ (k) - - (1 + V ) c« t (k) 
i=l ^ 

+ 53*[«-o; I (k)]^W)(k)-- C « t (k) 



53 % - w ,(k)] { ff ("'0(k) - ^4%(k)} = 
i=l 



(C6b) 
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The solution of these equations gives the following ex- 
pressions for the GF 



CT^)(k) 



^— ' uj — uJiik) ± id 



(C7a) 



G^(k, w ) = r(k) 

+ E a(rM) ( k ) 



i 



n 



uj + i 5 uj — i J 
0[w,(k)] 



a; — (k) + i S uj — ui (k) — i 5 _ 

(C7b) 



C^ t (k,«) =27r5(«)r(k) 

+ 2tt E S[uj-uji(k)]e[u>i(k)]a^(k) (C7c) 

C^(k,w) = 27r5(w)r(k) 

+ 2ttt7 E 5[w-^(k)]6»[-^(k)]CT("^)(k) (C7d) 

ig.A(k) 

It is fairly easy to check that these expressions (|C7ll cor- 
respond to limit T — > (/3 — > oo) of the expressions 
l|2TT3j) . 

APPENDIX D: USEFUL RELATIONS 

We note the dispersion relations 



G% A (k,oj) 

IT 



duj'- 



G ( g A (k,oj') } (Dla) 



5ft 



7T 



do/ 1 



uj — tu' 1 — 77 e" 



7- s 



(Dlb) 



This latter relation is valid for causal fermionic GF [i.e., 
for r) = 1] only when T(k) = 0. 

For the retarded and advanced GF, which are analyti- 
cal functions satisfying the standard Kramers-Kronig re- 
lations (jDlaj) . we can establish a spectral representation 



Gg ) A (k,w) = Jduj 



, P (,,) (k,c') 



— uj' ± i(5 

where we introduced the spectral function 

n 

p^(k,uj)=yS[uj^uj l (k)]a^ l \k) 



1=1 



=F 



7T L 



(D2) 



(D3) 



A spectral representation for the causal GF can be es- 
tablished in the following form 



rj e" 



G^(k )W ) = ffi/fM ! - A . r 

(D4) 

This latter relation is valid for causal bosonic GF [i.e., 
for rj = —1] only when T(k) = 0. 
We also note the sum rule 



/3w 



/n 
duj P M (k, u) = E ^ ( "' (k) = I {r,) (k) 

This is a particular case of the general sum rule 

/n 
duj uj p p M (k, w) = wfOO ct( "' ( k ) 

= M^' p )(k) = e p J 

where M^ ,p '(k) are the spectral moments 



(D5) 



(D6) 



M ( "' p) (k) =T 



d" 



(D7) 



and the last equality in Eq. (|D6|I holds only when Eq. Ijl.21 
) also does. 

Finally, by exploiting the independence of c"' (k) on rj 
[see Ea. H2.10c[) ]. we have 

a(- 1 ^(k)=tanh^^ ( T( 1 ^(k) V/ G B(k) (D8) 

In absence of symmetry breaking, Eqs. (|2.12|l and JD8J, 
together with Eq. (|2.8|) . give 



i#>(k) =E^(k)tanh^^Q r /(k)7«(k) 



(D9) 



The independence of G(k, uj) on 77 [see Eq. <|2. 141 )] gives 



G%\k,oj) 



G^"{k,uj) 
In terms of spectral densities 

ff (i.0(k) ( r(- 1 -0(k) 



(D10) 



E %"^(k)] 
Je.4(k) 



J _|_ e -/3 u)j (k) 1 _ e -/3 (k) 



= 
(Dll) 
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